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4.6 Dispersion of supension in a steady shear flow 

Similar to heat, a passive solvent is also transported by convection and diffusion. When the 
concentration of the solvent is low, the flow affects the transport but the solvent does not 
affect the flow. We explain a systematic theory of enhanced diffusion by shear- a phenomenon 
hrst analysized by G.I. Taylor, hence is called the Taylor dispersion. 

4.6.1 The convective diffusion equation for a passive suspension 

The spreading of a passive solute in a flow is governed by the convective diffusion equation 



where C is the volume concentration of the solute, and D the diffusivity. The derivation is 
as follows. 

Consider a material volume V. Conservation of the solute mass requires: 



+ V-(qC)= J DV 2 C 



(4.6.1) 




(4.6.2) 



The left hand side can be written as 




(4.6.3) 



Assuming Fick's law 
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(4.6.4) 



we get 




(4.6.5) 



Since V is arbitrary, (4.6.1) follows. 
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4.6.2 Heuristic picture of shear-enhanced diffusion and scales es- 
timates 

If a streak of die is put initially straight aross a channel of fluid in steady uniform flow, after 
some time the die streak will become curved due to the transverse variation of velocity, and 
also becomes thicker due to diffusion in all directions. See figure 4.6.1. After sufficiently 
long time lateral diffusion will be complete and the die distribution is nearly uniform across 
the channel. Convective diffusion will continue as a one-dimensional process. We wish to 
predict the effective diffusivity in the last phase. Consider dispersion in a one- dimensional 
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Figure 4.6.1: Initial stage of convective diffusion of a die streak. From Fischer et al. 

shear flow in a long circular pipe resulting from a pressure gradient which is constant in 
space but has steady and oscillatory parts. The velocity profile is known, 



u = U s (r), < r < a. 
The concentration of a solvent is governed by 

dC_ d(uC) _ lc?C_ ld_ ( dC\ 
dt dx \ dx 2 r dr \ dr I 



(4.6.6) 



(4.6.7) 



There are three time scales : diffusion time across the pipe radius a, convection time across 
L, and diffusion time across L. Their ratios are : 
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(4.6.8) 



The first time scale for lateral diffusion is very small. Let us focus attention to processes 
in the longitudinal direction long after lateral diffusion is complete. Hence we choose the 
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following normalization, 



x = Lx', r = or', u = U a u', t = ^-t', (4.6.9) 



where U Q is the scale of U. Equation (4.6.7) is normalized to 

U aa (dC d(u'C')\ a 2 d 2 C 1 d ( dC\ 



D L \dt' dx' I L 2 dx' 2 r dr \ dr 

Let the Peclet number Pe = U Q a/D be of order unity and a/L — e <C 1. Equation 
(4.6.10) becomes 

„ (dC d{u'C')\ 2 d 2 C 1 3 ( ,dC'\ 
with the boundary conditons 

ciC" 

— = 0, r' = 0,l (4.6.12) 

with 

v! = U' s (4.6.13) 
For brevity we drop the primes from now on. 

4.6.3 Multiple scale analysis-homogenization 

Ignoring the early stage of lateral difrusion, there are two time scales : convection time across 
L, and diffusion time across L. Their ratios are : 

— : — = 1 : - (4.6.14) 

Therfore we introduce the multiple time coordinates 

t, h = et, (4.6.15) 
and consider C(x, r, t, ti) to depend on t, and ti as if they were independent. In particular, 

dC dC dCdh dC dC 

-m^^ + df 1 ^ = ^ + e dT 1 (46 ' 16) 

Introduce multiple scale expansions. 

C = C + eC 1 + e 2 C 2 + .... (4.6.17) 
with Cj = Ci(x, r, t, ti), for alH = 1, 2, 3 • • •. The perturbation problems are 



0(6°): 

= l*L ( r ^ 

r dr \ dr 

with the boundary conditions: 

dC ° 0, r = 0,l. 



0(e): 



dr 



. ,'dCo d(uC )\ 1 d ( dC 
Pe — — + 



with: 



dt dx j r dr \ dr 
dC\ 



dr 



= 0, r = 0,l. 



0(e 2 



with 



dti dt dx J dx 2 r dr 



dC 2 

2 =0, r = 0,l. 



<9r 

The solution at 0(e°) is 

C = C (x,t 1 ) 

At 0(e), let the known velocity be 

u = U s (r) 

we get 



dt dx I r dr \ dr 

with 

dC x 

^ = 0, r = 0,l 
dr 

Let us now integrate(or average ) (4.6.26)across the pipe, and get 

dt + {Us) dx 

where angle brackets denote averaging over the cross section. 

1 r 1 
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Now subtract Pe(4.6.28) from (4.6.26) 



2nrh dr 

o 



p - dc id ( dc, 

PeU q — — = -— r 



dx r dr \ dr 
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where 

U = U s (r) - (U s ) (4.6.30) 

is the velocity nonuniformity 

Now C\ is governed by a linear equation, we can assume the solution to be proportional 
to the forcing, i.e., 

BC 

C 1 = Pe^-B s (r) (4.6.31) 

then 

Id/ d,B s \ ~ , . , 

with the boundary conditions 

^ = 0, r = 0,l. (4.6.33) 
dr 

and 

After solving for B s we go to 0(e 2 ), i.e., (4.6.22) . Note that , in view of (4.6.31), 



Pe— -i = Pe 2 «C/ s > + E/ S )£ s (r) 



ax ax" 1 

hence, 

„ /9Co dCA n2 ,, TT , ~ NrWN a 2 C <9 2 C 1 d / 0C 2 , 



which is a linear PDE for C*2. From (4.6.31) and (4.6.28) we find 

= -Pe 2 ^ 
at ax 2 



Pe ^l = _ Pe ^ (t/ , )B<(r . ) (4.6.35) 



P^ + W^-^ + i£fr^] (4,,0, 



It follows that 

Pp . , _ . , 

dt\ dx 2 dx 2 r dr \ dr 

with 

<K 2 r = 0, 1 (4.6.37) 



<9r 

Averaging (4.6.36) across the pipe, we get 
with 

£ = 1 - Pe 2 (U s B s ) (4.6.39) 

which is the effective diffusion coefficient. The first part is of molecular origin; the second 
part (the dispersion coefficient) is due to fluid shear. 
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Finally we add e (4.6.28) and (4.6.38) to get: 



This describes the convective diffusion of the area averaged concentration, which is certainly 
of practical value. 

After the perturbation analysis is complete, there is no need to use multiple scales; we 
may now write 

dC , x dC d 2 C 

-l + Pe(Us) -l = eE -l (4.6.41) 

In dimensionless form, it reads, 

This equation governs the convective diffusion of the cross-esctional average, after the 
initial transient is smoothed out. 

4.6.4 The dispersion coefficient for steady flow 

We demonstrate the details for a pure steady flow. 
The dimensional velocity profile is 

U s (r) = U (l- "pj (4.6.43) 
with U begin the maximum at the centerline. It is easy to show that 

(U) = ^ (4.6.44) 
We shall use U as the velocity scale for normalization, then in dimensionless form 

u(r) = 1 - r 2 , (u) = I (4.6.45) 



and the nonuniform part is 
The equation for B s is 



u = i - r 2 (4.6.46) 
1 d ( dB s \ 1 2 

-— r— - = - - r 2 4.6.47 

rdr\ dr 2 v ' 



which can be easily integrated to give 
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To detemine the integration constant we impose the condition that 

{B s ) =2 f 1 rB s dr = 
Jo 



This implies 



hence 



°^' + 24 



Using (4.6.46) and (4.6.51), we get 
<^,)=2^'r(i 

Thus 



— r 



24 8 



1 1 

h - I r" 

24 8 I 4 



(4.6.49) 

(4.6.50) 
(4.6.51) 

dr = -^ (4.6.52) 



9C ° + PeW dC ° 



In physical variables 



BC ° + { u,) dc ° 



at 



dx 



dx 



D 



e 1 + 



Pe 2 \ 9 2 Cn 



192/ 9x 2 



1 + 



192 V D 



d 2 C 
dx 2 



Denote by K the shear-enhanced dispersion coifficient, 



K = 



1 (U aY 
192 D 



(4.6.53) 



(4.6.54) 



(4.6.55) 



If the Peclet number is not small, the enhanced diffusivity (dispersivity) can be greater than 
the molecular diffusivity. Since the dispersivity is inversely proportional to D, its importance 
becomes greater for smaller molecular diffusisivity. 

Example: For salt in water, D = 10~ 5 cm 2 /s. If U Q = lem/s, and a = 0.2cm, then 
K = 22cm 2 /s, far exceeding the molecular diffusivity. 

Homework: Show that the shear-enhanced dispersion coefficient E for the uniform flow in 
a wide and inclined channel of depth h is (8/945)U 2 h 2 / D. 

Homework: Find the dispersion coefficient for sediments in an open channel flow. The 
sediments are uniform with fall velocity W Q . 



